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An atomic gas subject to a commensurate periodic potential generated by an optical lattice under- 
goes a superfluid-Mott insulator transition. Confining a strongly interacting gas to one dimension 
generates an instability where an arbitrary weak potential is sufficient to pin the atoms into the 
Mott state; here, we derive the corresponding phase diagram. The commensurate pinned state may 
be detected via its finite excitation gap and the Bragg peaks in the static structure factor. 



Atomic gases are developing into the ultimately tun- 
able laboratory system allowing to study complex quan- 
tum phenomena JjJ. Recently, subjecting an atomic 
Bose-Einstein condensate to an optical lattice, Greiner 
et al. Q have succeeded in tuning the system through 
a quantum phase transition separating a superfluid (S) 
from a Mott insulating (MI) phase. This 3D bulk tran- 
sition involves weakly interacting bosons and is well un- 
derstood within the Bose-Hubbard description J3|, |J : the 
system turns insulating when the on-site interaction en- 
ergy U becomes of the order of the hopping energy J. 
This strong coupling transition is a result of quenching 
the kinetic energy by a strong lattice potential. Amaz- 
ingly, by confining the atomic gas to one dimension, the 
strong coupling limit can be reached without the opti- 
cal lattice: in ID, the ratio 7 between the interaction- 
and kinetic energies per particle scales inversely with the 
density n and thus it is the Zoii>-density limit which is in- 
teracting strongly (Tonks gas) A new instability then 
appears in the strongly interacting ID quantum gas: the 
superfluid groundstate in the homogeneous system turns 
insulating in the presence of an arbitrarily weak commen- 
surate optical lattice and the S-MI transition changes to 
a transition of the incommensurate-commensurate type. 
In this letter, we analyze this new instability and derive 
the complete phase diagram for the S-MI transition in 
the limits of both weakly and strongly interacting gases. 
Remarkably, this goal can be achieved by a mapping to 
two classic problems, the Bose-Hubbard model (intro- 
duced in Ref. ^J) and the sine-Gordon problem describing 
the weakly- and strongly interacting limits of the atomic 
gas, respectively. Below, we first summarize the main 
results providing us with the phase diagram shown in 
Fig. 1; we then proceed with a detailed analysis of the 
strongly interacting Bose gas subject to a weak optical 
potential, leading us to a proper understanding of the 
new instability. 

A weakly interacting atomic gas subject to an opti- 
cal lattice is well described by the Bose-Hubbard model, 
which starts from a tight-binding model and takes the 
interaction between bosons into account perturbatively; 
the hopping amplitude J(V) and the on-site interaction 
energy U(V,j) follow from the underlying parameters of 
the atomic gas, the dimensionless interaction parameter 
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FIG. 1: Left: Schematic phase diagram illustrating the super- 
fluid and Mott insulating phases versus parameters 7 (interac- 
tion), V (optical potential), and Q (commensuration) . Right: 
Critical amplitude V c versus interaction 1 /j for the commen- 
surate situation with Q — 0. Below I/70 an arbitrary weak 
potential V drives the superfluid into the pinned insulating 
state. The dashed line denotes the asymptotic behavior near 
the critical point 1/7 C as determined from the sine-Gordon 
model, while the dashed-dotted line derives from the Bose- 
Hubbard criterion U/J\ S -mi ~ 3.84. 

7 and the amplitude V of the optical lattice. The phase 
diagram of the Bose-Hubbard model is well known 
and involves insulating Mott-lobes embedded in a super- 
fluid phase. In 3D, the mean-field analysis for densities 
commensurate with the lattice provides the critical pa- 
rameter U/J\ S -Mi ~ 5.8 z, in good agreement with the 
experimental findings of Greiner et al. S (here, z de- 
notes the number of nearest neighbors). Going to ID, 
fluctuations become important and appreciably modify 
the mean-field result: numerical simulations |7) place 
the transition at the critical value UJ J| S -mi = 2C ~ 3.84. 
This result is easily transformed into the j-V phase di- 
agram of the weakly interacting atomic gas, once the re- 
lations J(V) and U(V, 7) to the Mott-Hubbard param- 
eters are known; combining J oc exp[— 2(V/E r ) 1 / 2 ] as 
derived from a WKB-calculation, with the on-site repul- 
sion U oc 7, the transition line V c {^) (see Fig. 1) then 
derives from the implicit equation 

AV/E r = In 2 [4V2wC (V/E r ) 1/2 / 7 ] . (1) 

Here, E r = H 2 k 2 /2m is the recoil energy with m the bo- 
son mass and k = 2tt/X the wave vector of the light gener- 
ating the optical lattice V(x) — Vsm 2 (kx). The dimcn- 
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sionless interaction strength 7 is defined via 7 = mg/h 2 n, 
with n the density and g the strength of the 5-function 
interaction potential; g is related to the 3D scattering 
length a s and the transverse confining frequency uj±_ via 
g = 2hu> ± a s Q. 

Increasing the interaction strength 7, the critical am- 
plitude V c of the optical lattice triggering the S-MI tran- 
sition decreases, see Fig. 1; the description of the atom 
gas in terms of the Bose-Hubbard model breaks down 
and we have to look for a new starting point. For a weak 
optical potential, a natural choice is the ID Bose gas 
with ^-function interaction, which resides in the strong 
coupling regime at small densities n [jsj; the presence of 
the optical lattice is taken into account perturbatively. 
The homogeneous ID Bose gas with (5-function interac- 
tion has been solved exactly by Lieb and Liniger (|] ; the 
corresponding low-energy physics is properly described 
in terms of a Luttinger liquid with a parameter if (7) de- 
rived from the exact solution, see below. Adding the op- 
tical potential oc V, we arrive at the sine-Gordon model; 
the critical value K(j,V) = 2(1 + V/4E r ) separating 
the Mott insulating phase from the superfluid one de- 
termines the phase line V c /E r = (7^ — 7~ 1 )/5.5, with 
the critical value j c w 3.5. Combining the results of 
the Bose-Hubbard- and sine-Gordon models we can com- 
plete the phase diagram for the commensurable situation 
as shown in Fig. 1: most remarkable is the appearance 
of a critical interaction strength j c above which an arbi- 
trary weak optical lattice is able to pin the system into 
a Mott insulator state. The presence of this instabil- 
ity is due to the closeness of the dilute ID Bose liquid 
to Wigner-crystallization. Tuning the system away from 
commensurability with Q = 2ir(n — 2/A) 7^ the Mott in- 
sulator survives up to a critical misfit Q C (V, 7), see Fig. 1. 
The corresponding physics is similar to that of the com- 
mensurate-incommensurate transition of adsorbates on a 
periodic substrate as studied by Pokrovsky and Talapov 
& 

For a qualitative understanding of the transition in the 
strongly interacting ID Bose gas, it is useful to consider 
the limit 7 > 1, where the behavior is essentially that 
of an ideal Fermi gas pj]] . A weak periodic potential 
V(x) — Vsva^ikx), with a lattice constant such that 
an integer number i — 1,2,... of particles will fit into 
one unit cell, gives rise to a single-particle band struc- 
ture in which the i lowest bands are completely filled. 
The ground state for noninteracting Fermions is then 
a trivial band insulator, separated from excited states 
by an energy gap which scales like (V/E r y in the limit 
V -C E r . Similar to the Mott phase in the Bose-Hubbard 
model, the insulating state has a fixed integer density, 
commensurate with the lattice. It remains locked in a 
finite regime of the chemical potential, characterizing an 
incompressible state |3| . Clearly, for a weak periodic po- 
tential, the lowest energy gap with i = 1 is much larger 
than the higher-order ones. In the following, we will 



thus confine ourselves to studying the commensurate- 
incommensurate transition near integer filling i = 1, 
where the commensurate phase has maximal stability. 

For a quantitative theoretical analysis, it is convenient 
to use Haldane's description of the interacting ID Bose 
gas in terms of its long wave length density oscillations 
p2[ . Introducing the two fields cj>(x) and 9(x) describing 
phase and density fluctuations, the Hamiltonian of the 
homogeneous gas (without longitudinal confining trap 
and periodic potential) is a quadratic form involving ki- 
netic and interaction energies, 
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Here, vj = irhn/m is the analog of the bare 'Fermi' ve- 
locity at given average density n while vn = dn^i/irh is 
determined by the inverse compressibility, giving a sound 
velocity v s = y/vjVN consistent with the standard ther- 
modynamic relation mv 2 — nd n /i. For short-range repul- 
sive interactions, the dimensionless ratio K — (3 2 /An — 
vj/v s is larger than one, approaching 00 in the noninter- 
acting case and unity in the hard core limit fll2]| . In the 
idealized model with (5-function interactions, the exact 
solution by Lieb and Liniger shows that if is a mono- 
tonically decreasing function of the ratio 7 = mg/H 2 n 
between the interaction and kinetic energies; the limit- 
ing behavior for small values of 7 

if(7^0)= 7 r[ 7 -(l/2^)7 3 / 2 ]- 1 / 2 (3) 

follows from the Bogoliubov approximation in ID. Sur- 
prisingly, this result remains quantitatively correct for 7 
values up to 10 |J. At large 7 > 10, the asymptotic 
behavior is if (7 — > 00) = (1 + 2/^) 2 . 

For our subsequent analysis, it is convenient to intro- 
duce the conjugate fields <d(x) = (2/(3) [6(x) — irnx] and 
n(x) = —h(/3/2Tr)d x <f>(x) obeying the commutation rela- 
tion [0(x),n(x')] = ihS(x — x'). Using these fields, 

n„ = ^ / dx \(u/hf + (d x e) 2 } (4) 



takes the form of the Hamiltonian for a ID harmonic 
string with a linear spectrum to — v s q describing the 
long wave length density modulations of the interacting 
Bose gas. The assumption of a linear spectrum inherent 
in the simple form (||) is valid only below a momentum 
cutoff 1/a ~ to H ; the choice of the length scale a fixes 
the energy scale of Hq- The Boson density operator n(x) 
is related to the field 9(x) by 111 
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1 + 2 cos I — - — h Innx 
1=1 ^ 



(5) 



the last factor accounts for the discrete nature of the par- 
ticles. Adding an external periodic potential with ampli- 
tude V/2 and period A/2 gives rise to the perturbation 



V f Attx 
H v = — / dx n(x) cos — ^— . 
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As noted already by Haldane, insertion of the Fourier ex- 
pansion (|^) generates terms of the type appearing in the 
quantum (l+l)-dimensional sine-Gordon theory ]l3|, [li) . 
Close to commensurability the dominant term arising 
from the lowest harmonic in (|5]) has the conventional 
sine- Gordon form |0] 



V n f 

H w = — — j dx cos [/39 + Qa 
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with coupling parameter f3 = 2(nK) 1 ^ 2 and a twist 
Q = 2n(n — 2/ A). The strength of the nonlinear cos (3Q- 
perturbation is conveniently expressed through the di- 
mensionless parameter u = Tra 2 nV/2hv s which naturally 
involves the cutoff parameter a [l5| . The twist Q vanishes 
at commensurability; away from commensurability the fi- 
nite twist Q acts as a chemical potential for excitations 
and is preferably incorporated into the free Hamiltonian 
([|) via the replacement d x @ — > d x ® — Q//3. 

At fixed potential V, the quantum sine-Gordon model 
describes the competition between the preferred average 
inter-particle distance at given density due to the repul- 
sive interaction and the period imposed by the external 
potential. A perturbative calculation (see Ref. [l4| for a 
review) tells that for 1 /4tt = K > 2 a weak periodic po- 
tential is unable to pin the density; hence for K > K c = 2 
(7 < 7c ~ 3.5) the ground state remains gapless and su- 
perfluid in the presence of a small-amplitude lattice po- 
tential. In the strong coupling regime K < 2, however, 
the atoms are locked even to a weak lattice, as long as 
the twist Q is less than a critical value Q c . Beyond that, 
there is a finite density of 'solitons' (or domain walls for 
adsorbates on a periodic substrate, see Ref. |l^), which 
interpolate between minima of the external potential, re- 
lieving the frustration present at incommensurate densi- 
ties Q 7^ 0. The 'solitons' behave like relativistic parti- 
cles with energy E q — Tiv s \/ q 2 + M 2 and reestablish the 
superfluid response. The 'mass' M determines the exci- 
tation gap in the Mott insulating state, which translates 
into a jump A/i in the chemical potential at the commen- 
surate density: given that an additional/missing atom in- 
volves K solitonic excitations with energy E q= Q = hv s M 
|L6| one obtains 



An = 



2nfi 2 n 



in 



M: 
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furthermore, this mass is also simply related to the crit- 
ical twist via Q c = 2K 2 M. The precise numerical value 
of M depends on the high momentum cutoff 1/a via the 
dimensionless amplitude u of the lattice potential. The 
free-fermion limit K = 1 fixes this cutoff at I/a = nn, 
resulting in the simple form u = V/AE r ; we ignore small 
corrections arising due to a possible modification in the 
cutoff away from K = 1. The dependence of the mass 
M on u can be obtained from a recent nonperturba- 
tive renormalization group analysis of the quantum sine- 
Gordon model by Kehrein Jl5] ; for small values V <C E r 



and K away from K c 
takes the form 

An = 2E r 



the gap in the chemical potential 



V 



1 l/(2-K) 



(9) 
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In the limit K — > 1, the size of the gap approaches V/2, in 
agreement with the above fermionic picture of the Tonks 
gas limit, where the appearance of an insulating ground 
state, even for a weak periodic potential, is due to the 
opening of a single-particle band gap at the Fermi energy. 
In the practically accessible regime of large but finite 
7 > 7 C the gap depends on V in the more complicated 



manner as given by 
approaches K c = 2 fl5|. 



and vanishes exponentially as K 
see Fig. 2(a). Similarly, the den- 
sity range n — 2/X — ±Q c /2tt over which the ground state 
remains locked approaches zero as K — » K c . Finally, for 
K > 2, the dependence of the critical interaction param- 
eter K c (or 7 C ) on the lattice amplitude V follows easily 
from the Kosterlitz-Thouless nature of the scaling flow 
near K c = 2: to lowest order in u, K c (u) = 2(1 + u). 
Combining this result with (j^) it is straightforward to 
determine the line ^(7) separating the gapped insulat- 
ing regime from the superfluid at small but finite values 
of u, see Fig. 1. 



n 










FIG. 2: (a) Size of the gap Afi versus interaction strength 7 
for a fixed amplitude V = E r /2. For 7 — > 00 the gap assumes 
the free fermion limit V/2, while it vanishes exponentialy for 
7 — * 7c. (b) Fraction of atoms in the Mott insulating phase. 
The inset shows the density distribution n(x) with the Mott 
phase characterized by a locked commensurate density in the 
trap center, surrounded by a superfluid region. 

In order to analyze the consequences of the commensu- 
rate-incommensurate transition for cold atoms in a trap, 
we consider a ID Bose gas in the presence of a weak longi- 
tudinal confining potential V(x) = mto 2 x 2 /2. Provided 
the associated oscillator length I — (h/muj) 1 / 2 is much 
larger than the inter-particle distance, the density profile 
in this inhomogeneous situation may be obtained from 
the Thomas-Fermi approximation fl7| 

(j, [n(x)] + V(x) = /i[n(0)] (10) 

where /x[n] is the chemical potential of the homogeneous 
system. The central density n(0) and the associated ra- 

1 /2 

dius R = (2/i[?i(0)]/mcj 2 ) of the cloud is obtained 
from the normalization condition 



N = 



[ dx n{x) = 2R [ 
J-r Jo 



(0) 



dn\ / 1 — 



n[n(o)] 



(11) 



4 



with N the particle number, provided the relation /z[n] 
is known explicitly. 

In the absence of an optical lattice the density profiles 
are smooth functions of the coupling 7 |l7j . In the limit 
7> 1, [i[n] — * /j F [n] = (Tntn) 2 /2m approaches the Fermi 
energy of an ideal Fermi gas with density n, resulting 
in a profile n(x) = (2N/irR)y/l — (x/R) 2 with radius 
R = (2Ny/ 2 l. Adding now a periodic potential which is 
nearly commensurate with the density in the trap center 
opens a gap in the chemical potential. The density 
profile will then exhibit a flat, incompressible regime in 
the trap center. In order to find the detailed shape of the 
density profile, we approximate the sine-Gordon model in 
its strong coupling, gapped phase by a gas of noninter- 
acting relativistic Fermions of mass M; in the relevant 
regime 1 < K < 2, this approximation is known to work 
extremely well . The chemical potential as a function 
of density can then be calculated explicitly, 

r , Au / AKE r 8n\ , x 

mW«M»] + ^/(-^-J (12) 

with 5n = n—n c (n c = 2/ A). The dimensionless function 

f(z) = ±{l + z 2 ) 1/2 -z (13) 

is discontinous at z = and incorporates, via (p"2|), the 
jump Afi in the chemical potential at n c associated with 
the incompressible commensurate state. Using this ap- 
proximation, the density profile and the fraction of parti- 
cles participating in the commensurate phase derive from 
an integration of ( |TT| ) with (jl^), see Fig. |(b). Knowl- 
edge of the locked fraction AN/N plays an important 
role in the experimental detection of a commensurate 
Mott phase. This can be achieved by measuring the ex- 
citation gap through a phase gradient method as done 
previously for the Bose-Hubbard transition || . Alterna- 
tively, it should be possible to directly observe the in- 
crease in the long-range translational order in the Mott 
phase via Bragg diffraction []l8[ in either case the 
fraction AN/N determines the experimentally available 
signal. The latter can be further enhanced by generating 
an array of parallel 'atom wires' with the help of a strong 
2D optical transverse lattice. Using numbers similar to 
those in the recent experiment by Greiner et al. , it is 
possible to generate several thousand parallel ID wires 
with a transverse confining frequency v± = 20 kHz. A 
longitudinal harmonic trap with frequency v = 40 Hz 
then encloses N w 50 atoms in each ID wire; the as- 
sociated central density in the absence of a longitudinal 
periodic potential is n(0) = 2 fivaT 1 for 7 » 1, com- 
mensurate with the lattice constant A/2 « 0.5 /zm of a 
typical optical lattice j| . A weak periodic potential will 
then lead to an incompressible Mott state in the center 
of the cloud, provided the parameter 7 = 2a s /n(0)l 2 L is 



larger than the critical value 7 C « 3.5. For 87 Rb with 
a scattering length a s w 5 nm, the resulting 7 is equal 
to one, i.e., not quite in the required range. As noted 
already by Petrov et al. 0, however, larger and in par- 
ticular tunable values of 7 may be realized by changing 
a s via a Feshbach resonance as present, e.g., in 85 Rb. 

In conclusion we have shown that a commensurate 
Mott state can be realized in dilute ID BEC's already 
with an arbitrary weak lattice potential, provided that 
the ratio 7 between the interaction and kinetic energies 
is larger than a critical value 7 C w 3.5. This instability 
provides a new and experimentally accessible tool for the 
quantitative characterization of ID atomic gases in the 
strongly correlated 'Tonks gas' limit. Also, the obser- 
vation of a Mott state in a regime where the atoms are 
not confined to discrete lattice sites would give direct evi- 
dence for the granularity of matter in strongly interacting 
dilute gases. 
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